The low-density spin susceptibility and effective mass of mobile electrons 

in Si inversion layers 
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We studied the Shubnikov-de Haas (SdH) oscillations in high-mobility Si-MOS samples over a 
wide range of carrier densities n ~ (1 — 50) x 10 11 cm -2 , which includes the vicinity of the apparent 
metal-insulator transition in two dimensions (2D MIT). Using a novel technique of measuring the 
SdH oscillations in superimposed and independently controlled parallel and perpendicular magnetic 
fields, we determined the spin susceptibility \* \ the effective mass m* , and the (/'-factor for mobile 
electrons. These quantities increase gradually with decreasing density; near the 2D MIT, we observed 
enhancement of \* by a factor of ~ 4.7. 
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Many two-dimensional (2D) systems exhibit an appar- 
ent metal-insulator transition (MIT) at low temperatures 
as the electron density n is decreased below a critical den- 
sity n c (for reviews see, e.g., Refs. [Q— g| ) . The phenom- 
ena of the MIT and 'metallic' conductivity in 2D attract 
a great deal of interest, because it addresses a funda- 
mental problem of the ground state of strongly correlated 
electron systems. The strength of electron-electron inter- 
actions is characterized by the ratio r s of the Coulomb 
interaction energy to the Fermi energy ef- The 2D MIT 
is observed in Si MOSFETs at n c ~ 1 X 10 n cm~ 2 , 
which corresponds to r s ~ 8 (for 2D electrons in (lOO)-Si, 
r s = 2.63v/l0 12 cm- 2 /n). 

In the theory of electron liquid, the electron effec- 
tive mass to*, the g*-factor, and the spin susceptibility 
X* oc g*m* are renormalized depending on r s [Q. Though 
the quantitative theoretical results vary consider- 

ably, all of them suggest enhancement of x*, to* and g* 
with r s . Earlier experiments d[ll]] have shown growth of 
to* and g*m* at relatively small r s values, pointing to a 
ferromagnetic type of interactions in the explored range 
1 < r s < 6.5. Potentially, strong interactions might drive 
an electron system towards ferromagnetic instability [Q . 
Moreover, it has been suggested that the 'metallic' be- 
havior in 2D is accompanied by a tendency to a ferro- 
magnetic instability [Q. Thus, in relation to the still 
open question of the origin of the 2D MIT, direct mea- 
surements of these quantities in the dilute regime near 
the 2D MIT are crucial. 

In this Letter, we report the direct measurements of 
X*, m*, and hence g* over a wide range of carrier den- 
sities (1 < r s < 8.4), which extends for the first time 
down to and across the 2D MIT. The data were obtained 
by a novel technique of measuring the interference pat- 
tern of Shubnikov-de Haas (SdH) oscillations in crossed 
magnetic fields. The conventional technique of measur- 
ing g*m* in tilted magnetic fields [p|JTc|] is not applicable 



when the Zeeman energy is greater than the cyclotron 
energy [Q. The crossed field technique removes this re- 
striction and allows us to extend measurements over the 
wider range of r s . We find that for small r s , the g*m* 
values increase slowly in agreement with the earlier data 
by Fang and Stiles || and Okamoto et al. For larger 
values of r s , g*m* grows faster than it might be extrap- 
olated from the earlier data. At our highest value of 
r s ~ 8.4, the measured g*m* is greater by a factor of 4.7 
than that at low r s . 

Our resistivity measurements were performed by ac 
(13Hz) technique at the bath temperatures 0.05 — 1.6 K 
on six (100) Si-MOS samples selected from 4 wafers: 
Sil2 (peak mobility fx ~ 3.4m 2 /Vs at T = 0.3 K), Si22 
(p ~ 3.2m 2 /Vs), Si57 (p ~ 2.4m 2 /Vs), Si6-14/5, Si6- 
14/10, and Si6-14/18 (p ~ 2.4m 2 /Vs for the latter three 
samples) jt4|]. The gate oxide thickness was 190 ± 20 nm 
for all the samples; the 2D channel was oriented along 
[011] for sample Si3-10, for all other samples - along [010]. 
The crossed magnetic field system consists of two mag- 
nets, whose fields can be varied independently. A split 
coil, producing the field B± < 1.5T normal to the plane 
of the 2D layer, is positioned inside the main solenoid, 
which creates the in-plane field B\\ up to 8T. The elec- 
tron density was determined from the period of SdH os- 
cillations. 

Typical traces of the longitudinal resistivity p xx as a 
function of B± are shown in Fig. 1. Due to the high 
electron mobility, oscillations were detectable down to 
0.25 T and a large number of oscillations enabled us to 
extract the fitting parameters g*m* and to* with a high 
accuracy. The oscillatory component 8p xx was obtained 
by subtracting the monotonic 'background' magnetore- 
sistance (MR) from the p xx (B±) dependence (the back- 
ground MR is more pronounced for lower n, compare 
Figs. 1 and 3). 

The theoretical expression for the oscillatory compo- 
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nent of the magnetoresistance is as follows pi 
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Here po — p xx (B± = 0), u> c = eB±/m*m e is the cy- 
clotron frequency, m* is the dimensionless effective mass, 
m e is the free electron mass, and Tjj is the Dingle tem- 
perature. We take the valley degeneracy g v = 2 in 
Eqs. (1), (2) and throughout the paper. The Zeeman 
term Z s = cos[7rs?i7r(n| — n^)/ (eB±)] reduces to a field- 
independent constant for £?|| = 0. 

Application of Bu induces beating of SdH oscillations, 
which are observed as a function of B±. The beat fre- 
quency is proportional to the spin polarization of the 



interacting 2D electron system |16 
X*B tot 
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where n-\ (nj.) stands for the density of spin-up (spin- 
down) electrons, gt — 2 is the bare g-factor for Si, and 

B tot = <Jb\ +-Bjj. Equations (2), (3) imply that the 
spin polarization of the electron system is linear in B to t 
(relevance of this assumption to our experiment will be 
justified below). 
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FIG. 1. Shubnikov-de Haas oscillations for 
n = 10.6 x 10 n cm- 2 (sample Si6-14/10) at T = 0.35K 
and two values of Bu. 

In the experiment, we observed a well pronounced 
beating pattern at a non-zero Bu (Figs. 1 and 2), in 
agreement with Eq. (1). The phase of SdH oscillations 
remains the same between the adjacent beating nodes, 
and changes by ir through the node. The positions of 



the nodes are defined solely by g*m*; observation of sev- 
eral nodes enables us to determine this quantity with a 
high accuracy. The positions of the nodes on the B± 
axis (and, thus, the value of g*m*) are T-independent at 
T < IK within the experimental accuracy ~ 2%. We 
have observed a non-monotonic dependence of g*m* on 
which will be discussed elsewhere p7j . This depen- 
dence is more pronounced near the 2D MIT, where g*m* 
varies with B\\ by ~ 15%. The data discussed below have 
been obtained in the fields sufficiently low to ignore the 
g*m*(B) dependence. For example, to ensure the linear 
regime, at high densities (n ~ 10 12 cm~ 2 ) we have used 
B\\ ~ 0.8 T , which corresponded to the spin polarization 
P ~ 1%; at low densities (n ~ 10 n cm~ 2 ), the applied 
Sii ~ 0.1 - 0.2 T corresponded to P - 5 - 10%. 
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(a) m* = 0.22, T D = 0.55K, g* = 2.63 
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FIG. 2. Examples of data fitting with Eq. (1): (a) 
n = 10.6 x 10 u cm^, T = 0.35 K, By = 4.5 T, P = 6%; 
(b) n = 2.2xlO n cm- 2 , T = 0.4 K, By = 3.36 T, P = 34%. 
The data for sample Si 6-14/10 are shown as the solid lines, 
the fits (with parameters shown) as the dashed lines. Both 
are normalized by Ai(B±). 

Comparison between the measured and calculated de- 
pendences Sp xx /po versus B±, both normalized by the 
amplitude of the first harmonic A\ , is shown in Fig. 2 for 
two carrier densities. The normalization assigns equal 
weights to all oscillations. We analyzed SdH oscilla- 
tions over the low-field range B± < IT; this limitation 
arises from the assumption in Eq. (1) that Hoj c < ej? 
and 5p xx /po <C 1. The latter condition also allows us 
to neglect the inter-level interaction which is known to 
enhance g* in stronger fields pSy . 

The amplitude of SdH oscillations at small B± can be 
significantly enhanced by applying Bu (see Fig. 3), which 
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is another advantage of the cross-field technique. Indeed, 
for low B± and n, the electron energy spectrum is com- 
plicated by crossing of levels corresponding to different 
spins/ valleys. By applying Sm, we can control the energy 
separation between the levels, and enhance the amplitude 
of \ow-B oscillations (see Fig. 3). We have verified that 
application of By (up to the spin polarization ~ 20%) 
does not affect the extracted m* values (within 10% ac- 
curacy), provided the sample remains in the 'metallic' 
regime (the insets to Fig. 3 show that the values of to* 
measured at Bn = and 3.36 T do coincide). 
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FIG. 3. Shubnikov-de Haas oscillations versus B± 
sample Si 6-14/10 (n = 2.2 x 10 n cm" 2 ) at T = 0.4; 0.5; 
0.6; 0.7; 0.8K: (a) B\\ = and (b) B\\ = 3.36 T. The insets 
show the temperature dependences of fitting parameters 
(T + T D )m*. 

Fitting of the data provides us with two combinations 
of parameters: g*m* and (T + Tjj)m*. The measured 
values of g*m*, as well as to* which are discussed be- 
low, were similar for different samples. Figure 4 a shows 
that for small r s , our g*m* values agree with the earlier 
data by Fang and Stiles j|] and Okamoto et al. |l0|. For 
r s > 6, g*m* increases with r s faster than it might be ex- 
pected from extrapolation of the earlier results |l(J . The 
MIT occurs at the critical r s -value ranging from r c = 7.9 
to 8.8 for different samples; in particular, r c ~ 8.23 for 
samples Si6-14/5,10,18 which have been studied down to 
the MIT. 

An interesting question is whether the measured de- 
pendence X*( r s) shown in Fig. 4a represents a critical 
behavior as might be expected if x* diverges at n c . Anal- 
ysis of the Bn -induced magnetoresistance led the authors 



of Refs. |2^,^3| to the conclusion that there is a ferro- 
magnetic instability at n = n c . By forcing the critical 
dependence x* K ( n / n x ~ l) -a to fit our data, we found 
that the correlated fitting parameters n x and a can be 
varied over a wide range, only weakly affecting the stan- 
dard deviation. More importantly, the period of SdH 
oscillations shows that the electron states remain double 
spin degenerate across the 2D MIT down to at least n 
— 0.98x 10 11 cnr 2 . The details of this analysis, which 
does not support the conclusion on spontaneous spin po- 
larization and divergency of x* & t n — n c are presented 
elsewhere lO. 
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FIG. 4. Parameters g*m* , m* , and g* for different sam- 
ples as a function of r s (dots). The solid line in Fig. 4a 
shows the data by Okamoto et al. |l(]]. The solid and 
open dots in Figs. 4 b and 4 c correspond to two differ- 
ent methods of finding m* (see the text). The solid and 
dashed lines in Fig. 4b are polynomial fits for the two de- 
pendences m*(r s ). The values of g* shown in Fig. 4c 
were obtained by dividing the g*m* data by the smooth 
approximations of the experimental dependences m*(r 3 ) 
shown in Fig. 4 b. 

The second combination, (T+To)m*, controls the am- 
plitude of oscillations. In order to disentangle Tjj and 
to*, we analyzed the temperature dependence of oscilla- 
tions over the range T = 0.3 — 1.6K (for some samples 
0.4 — 0.8 K) [^9|. The conventional procedure of calculat- 
ing the effective mass for low r s values ($5), based on 
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the assumption that Try is T-independent, is illustrated 
by the insets in Fig. 3. In this small-r s range, our results 
are in a good agreement with the earlier data by Smith 
and Stiles 0, and by Pan et al. jy]]. The assumption 
Tjj ^ f{T), however, becomes dubious near the MIT, 
where the resistance varies significantly over the studied 
temperature range; in this case, the two parameters Tjj 
and to* become progressively more correlated. The open 
dots in Fig. 4 b were obtained by assuming that To is 
T-independent over the whole explored range of n: m* 
increases with r s , and the ratio m* /m& becomes ~ 2.5 at 
r s — 8 (rrib = 0.19 is the band mass). As another limit- 
ing case, one can attribute the change in R(T) solely to 
the temperature dependence of the short-range scatter- 
ing and request To to be proportional to R(T). In the 
latter case, the extracted dependence m*(r s ) is weaker 
(the solid dots in Fig. 4b). To reduce the uncertainty of 
to* at large r s , it is necessary to separate the effects of 
T-dependent scattering and 'smearing' of the Fermi dis- 
tribution in the dependence R(T); the adequate theory 
is currently unavailable. 

Our data shows that the combination (T + T^)m* 
is almost the same for electrons in both spin-up and 
spin-down subbands (e.g., for n — 3.76 x 10 11 cm -2 and 
B\\ = 2.15 T (P = 20%), the T D values for 'spin-up' and 
'spin-down' levels differ by < 3%). This is demonstrated 
by the observed almost 100% modulation of SdH oscil- 
lations (see, e.g., Fig. 2b). Thus, the carriers in the 
spin-up and spin-down subbands have nearly the same 
mobility; this imposes some constraints on theoretical 
models of electron transport in the 2D 'metallic' state. 

In conclusion, using a novel cross-field technique for 
measuring SdH oscillations, we performed direct mea- 
surements of the spin susceptibility, effective mass, and 
g— factor of conducting electrons over a wide range of car- 
rier densities. By studying the temperature dependence 
of SdH oscillations, we disentangled three unknown pa- 
rameters: g* , m*, and To, and obtained their values in 
the limit of small magnetic fields. We found that both 
g* to* and to* are almost independent of the tempera- 
ture over the range (0.2 — 1) K. At the 2D metal-insulator 
transition, we observed a finite value of the spin suscep- 
tibility, enhanced by a factor of ~ 4.7 with respect to the 
high-density limit. 
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